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Phases of matter are conventionally characterized by order parameters describ¬ 
ing the type and degree of order in a system. For example, crystals consist of 
spatially ordered arrays of atoms, an order that is lost as the crystal melts. Like¬ 
wise in ferromagnets, the magnetic moments of the constituent particles align 
only below the Curie temperature, Tq. These two examples reflect two classes of 
phase transitions: the melting of a crystal is a first-order phase transition (the 
crystalline order vanishes abruptly) and the onset of magnetism is a second- 
order phase transition (the magnetization increases continuously from zero as 
the temperature falls below Tq). Such magnetism is robust in systems with lo¬ 
calized magnetic particles, and yet rare in model itinerant systemd30 where the 
particles are free to move about. Here for the first time, we explore the itiner¬ 
ant magnetic phases present in a spin-1 spin-orbit coupled atomic Bose gad^^, 
in this system, itinerant ferromagnetic order is stabilized by the spin-orbit cou¬ 
pling, vanishing in its absence. We first located a second-order phase transition 
that continuously stiffens until, at a tricritical point, it transforms into a first- 
order transition (with observed width as small as h x 4 Hz). We then studied 
the long-lived metastable states associated with the first-order transition. These 
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measurements are all in agreement with theory. 


Most magnetic systems are composed of localized particles such as electron^, atomic 
nucleP, and ultracold atoms in optical latticeP^, each with a magnetic moment fi. By 
contrast, itinerant magnetisnP describes systems where the magnetic particles - here ultra¬ 
cold neutral atoms - can themselves move freely, and for which magnetism is generally weak. 
Our spin-orbit coupled Bose-Einstein condensates (BECs) constitute a magnetically ordered 
itinerant system in which the atoms’ kinetic energy explicitly drives a phase transition be¬ 
tween two different ordered phaseP. The coupling between spin and momentum afforded by 
spin-orbit coupling (SOC) provides a new route for stabilizing ferromagnetism in itinerant 
systems, both for fermionJ^ and here for boson^lH 

Phase transitions can generally be described in terms of a free energy G{Mz) - including 
the total internal energy along with thermodynamic contributions - that is minimized for 
an equilibrium system. Here the magnetization = {Sz)/h is an order parameter which 
changes abruptly as our system undergoes a phase transition, where S is the spin. Figure 
shows typical free energies: a hrst-order phase transition (top panel) occurs when the number 
of local minima in G{Mz) stays hxed, but the identity of the global minima changes; and a 
second-order phase transition (bottom panel) occurs when degenerate global minima merge 
or separate. These dehning features are true both for T > 0 thermal and T = 0 quantum 
phase transitions. Here we realize / = 1 BECs with SOC and study the magnetic order and 
associated quantum phase transitions present in their phase diagram. 
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For spin-1/2 systems (i.e, total angular momentum, / = 1/2) like electrons, ferromag¬ 
netic order can be represented in terms of a magnetization vector M = {S)/h. This is rooted 
in the fact that the three components of the spin operator S transform vectorially under 
rotation. More specihcally, any Hamiltonian describing a two level system may be expressed 
a,s H = hQo + the sum of a scalar (rank-0 tensor) and a vector (rank-1 tensor) contri¬ 

bution. The former, described by Hq gives an overall energy shift, and the latter takes the 
form of the linear Zeeman effect from an effective magnetic held proportional to fii. Going 
beyond this, fully representing a spin-1 (i.e., total angular momentum / = 1) Hamiltonian 
with angular momentum F rather than spin S, requires an additional hve-component rank-2 
tensor operator - the quadrupole tensor - and therefore there exist “magnetization” order 
parameters that are not simply associated with any spatial directionP^ISI 

Pioneering studies in GaAs quantum well£^E3 showed that material systems with equal 
contributions of Rashba and Dresselhaus SOG described by the term 2hkRkxFz/m, subject 
to a transverse magnetic held with Zeeman term can equivalently be described as a 

spatially periodic ehective magnetic held. Our spin-orbit coupled spin-1 atomic system^MUlEll 
have this form of SOG and can therefore be described by the magnetic Hamiltonian 

H = ^ + n,{x)-F + n,F!:!K (i) 

describing atoms with mass m and momentum Kk interacting with an ehective Zeeman 
magnetic held r2i(x)/r2i = cos{2kjix)ex — sm{2kfix)ey helically precessing in the e^-Sy plane 
with spatial period vr/fc^ set by the SOG strength; and an additional Zeeman-like tensor 
coupling with strength 02. Here, F^^/h = .F//h^ — 2/3 is an element of the quadrupole 
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tensor operator. 


The competing contribntions - from kinetic and magnetic ordering energies - make 
onrs an archetype system for stndying robnst itinerant magnetic order and nnderstanding 
the associated phase transitions, of which both hrst- and second- order are present (Fig. Hf)' 
For inhnitesimal f2i, the tensor held favors either: for f22 > 0, a polar BEC {mp = 0: 
nnmagnetized, = 0), or for f22 < 0, a ferromagnetic BEC {irtp = -|-1 or —1: magnetized, 
\Mz\ = 1); these phases are separated by a hrst-order phase transition at f22 = 0. In contrast, 
for large ffi the system forms a spin helix BEC (with local magnetization antiparallel to fii: 
nnmagnetized, = 0). This order increases the system’s kinetic energy, leading to the 
second-order phase transition shown in Fig. These two extreme limits continnonsly 
connect at the point the green star in Fig. [^, where the small-fli hrst-order 

phase transition gives way to the large-f2i second-order transition, and together these regions 
constitnte a cnrve of critical points {(f2f,fl^)}. 

As shown in Fig. we realized this sitnation by illnminating ^^Rb BECs in the 
total angnlar momentnm / = 1 gronnd state manifold with a pair of connter propagating 
“Raman” lasers that coherently conpled the manifold’s three mp states. As we hrst showecP 
nsing ehective / = 1/2 systems, this introdnces both a spin-orbit and a Zeeman term into the 
BEC’s Hamiltonian, eqnivalent to Eq. ([^. The Raman laser wavelength A dehnes the natnral 
nnits for onr system: the single-photon recoil energy Ep = h'^kp^/2m, and momentnm 
hkp = 27rh/A setting the SOC strength. Here the qnadratic Zeeman shift from a large 
bias magnetic held BqBz split the low-held degeneracy of the \mp = —1) <H- \mp = 0) and 
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\mp = Q) \mp = +1) transitions, and we independently Raman conpled these state- 
pairs with eqnal strength f2i. As described in the Methods Snmmary, we dynamically tnned 
the qnadrnpole tensor held strength f22 by simnltaneously adjnsting the Raman freqnency 
differences. Withont this techniqne, only the npper half-plane of the phase diagram (Fig. [^) 
wonld be accessible: containing only an nnmagnetized phase, therefore lacking any phase 
transitions. 


In each experiment, we first prepared BECs at a desired point in the phase dia¬ 
gram, possibly having crossed the phase transition dnring preparation; a combination of 
trap dynamic^25112§l, collisions, and evaporation!^ kept the system in or near (local) thermal 
eqnilibrinm. We then made magnetization measnrements directly from the spin resolved 


momentnm distribntion obtained nsing the time-of-ffight techniqnes described in Ref. 26 


Onr experiment first focnsed on thermodynamic phase transitions. We made vertical 
(horizontal) scans throngh the phase diagram by initializing the system in the nnmagnetized 
phase at a desired valne of fli (fl 2 ) with ^ 0 (fli < lOE/j), and then ramping f 22 
(fli) throngh the transition region. (As discnssed in the methods snmmary onr nominally 
horizontal f2i-scans followed slightly cnrved trajectories throngh the phase diagram, snch as 
the red dashed cnrve in Fig. [^). Following snch ramps, domains with both can rapidly 
form, and we therefore focns on the tensor magnetization M^z = {Fz^)/f^~\~ 2/3 which is 
sensitive to this local magnetic order. 


Using horizontal scans, we crossed throngh the second-order phase transition (f22 < ^^ 2 ) 
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where the free energy evolves continuously from having one minimum (with M^z = 0, for 
large hli) to having two degenerate minima (with Mzz > 0, for smaller hli). As shown 
in Fig. if. M^z continuously decreases, reaching its saturation value as hli exceeds hlf. 
Repeating this processes for < ^2 < 0, we found a sharp hrst-order transition. In each 
case, data is plotted along with theory with no adjustable parameters. Using data of this 
type for a range of 1^2 and htting to numeric solutions of Eq. ([^, we obtained the critical 
points plotted in Fig. if. Because horizontal cuts through the phase diagram are nearly 
tangent to the transition curve for small ^ 2 , this produced large uncertainties in hlf for the 
hrst-order phase transition. 

We studied the hrst-order phase transition with greater precision by ramping ^2 through 
the transition at hxed Ui (Fig. [^) and found near perfect agreement with theory. For all 
the experimentally measured critical points, see Fig. top, separating the unmagnetized 
and ferromagnetic phase, we also measured the corresponding transition widths. Dehned as 
the required interval for the curve to fall from 50% to 20% of its full range, this width A 
decreases sharply at marking the crossover between second- and hrst-order phase tran¬ 
sitions (see Fig. bottom). In these data, the width of the hrst-order transition becomes 
astonishingly narrow: as small as 0.0011 (3)^^ = h x 4(l)Hz at fli = 0.41(1). 

We observed that scans crossing the second-order transition typically required 50 ms 
to equilibrate, while for scans crossing the hrst-order transitions we allowed as long as 1.5 s 
for equilibration. Systems taken through a hrst-order phase transition can remain in long- 
lived metastable states. Here a metastable state with Mz = 0 persists in the ferromagnetic 
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phase, and a pair of metastable states with 7 ^ 0 persists in the nnmagnetized phase. 
We began onr study of this metastability by quenching through the hrst-order transition at 
fii = 0.74(8)ii^ij with differing rates from 0.5 to 0.2 Eji/s, as shown in Fig. We observed 
the transition width continuously decreases with decreasing ramp rate (inset to Fig. [^, 
consistent with slow relaxation from a metastable initial state. 

We explored the full regime of metastability by initializing BECs in each of the Imp = 
0 , ± 1 ) states, at hxed then rapidly ramping from zero to its hnal value while 

remaining adiabatic when changed rapidly (but still fast enough so that metastable 
states survive, < 200E^/s). For points near the hrst-order phase transition three metastable 
states exist (Fig. |^; near the second-order transition this count decreases, giving the two 
local minima which merge to a single minimum beyond the second-order transition. 

We experimentally identihed the number of metastable states by using and its 
higher moments, having started in each of the three \mF = 0 ,± 1 ) initial states. A small 
variance indicates the hnal states are clustered together - associated with a single global 
minimum in the free energy G{Mz) - and it increases when metastable or degenerate ground 
states are present. We distinguished systems with two degenerate magnetization states 
{Mz ~ ± 1 ) from those with three states by the same method, since when ~ ± 1 , the 
variance of \Mz\ is small, and it distinguishably increases as a third metastable state appears 
with Mz = 0. In this way we fully mapped the system’s metastable states in agreement with 
theory, as shown in Fig. 
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In conclusion, we accurately measured the two-parameter phase diagram of a spin-1 
BEC, containing a ferromagnetic phase and an unmagnetized phase, continuously connecting 
a polar spinor BEC to a spin-helix BEC. The ferromagnetic phase in this itinerant system 
is stabilized by SOC, and vanishes as the SOC strength hkn goes to zero. Our observation 
of controlled quench dynamics through a hrst-order phase transition opens the door for 
realizing Kibble-Zurek physic^l^ in this system, where the relevant parameters can be 
controlled at the individual Hz level. The quadrupole tensor held oc studied here is the 
q = 0 component of the rank-2 spherical tensor operator with q G {±2, ±1,0}. The 
physics of this system would be further enriched by the addition of the remaining four tensor 
helds: this is straightforward using a combination of radio-frequency and microwave heldd^. 
Furthermore, the hrst-order phase transition revealed in this work is further modihed by the 
addition of spin-dependent interactions: very weak in ®^Rb, giving ehects below our current 
ability to resolve. Our numerical studies show that a number of new symmetry-broken spinor 
phases emerge along the line that dehnes the hrst-order transition. 

Methods Summary 

System Preparation We created ^ 4 x 10® atoms ®^Rb BECs in the ground electronic 
state 1/ = 1) manifold^, conhned in the locally harmonic trapping potential formed at the 
intersection of two 1,064 nm laser beams propagating along and e^ giving trap frequencies 
of {ux,u}y,ujz)= (33(2),33(2), 145(5)) Hz. The quadratic contribution to the Bq = 
35.468(1) G magnetic held’s ^ h x 25 MHz Zeeman shift lifts the degeneracy between the 



1 / = 1 , mp = — 1 ) <H- \ f = l,mp = 0) and \f = 1 , mp = 0 ) -H- |/ = 1 , mp = 1 ) transitions, by 
e = h X 90.417(1) kHz. We denote the energy differences between these states as hSu-i^, 
h6uo^pi, and 

Frequency selective Raman coupling We Raman-coupled the three mp states using a 
pair of A = 790.024(5) nm laser beams counter-propagating along e^,. The beam traveling 
along +6^ had frequency components and ojXii while the beam traveling along —e^. 
contained the single frequency u:~. These frequencies were chosen such that the differences 
(5a;_i,o ~ W — and 5a;o,+i ~ W — independently Raman coupled the \mp = —1, 0) 
and \mp = 0,-|-1) state-pairs, respectively. Furthermore we selected 2a;“ — (cjI'^i Tcj+i) = 
such that, after making the rotating wave approximation (RWA) the \mp = ±1) 
states were energetically shifted by the same = [((5ti;_i,o — <^i^o,+i) + (i^-i — 
energy from \mp = 0), thereby yielding the frequency-tuned tensor energy shift depicted 
in Fig. [^. In addition, 112 in Eq. ([^ differs from by a small shift oc Hi/e resulting 
from off-resonant coupling to transitions detuned by 2 e, which we computed directly using 
Floquet theory [see Eq. ( 6 ) in Methods]. This tensor contribution VL 2 F^X’ to the Hamiltonian 
might equivalently be introduced by a quadratic Zeeman shift alone, giving H 2 oc Rg > 0. 
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Figure 1 | Experimental system a Schematic and level diagram. The Imp = - 1 ) -h- Imp = 
0 ) and \mp = 0) -H- \mp = + 1 ) transitions of the / = 1 ground state manifold of ^'^Rb were 
independently Raman coupled, giving experimental control of fli and ^ 12 - b Phase diagram. The 
ferromagnetic order parameter \M^\ is plotted against ^^2 and rii. The solid (dashed) red curve 
denotes the first-order (second-order) transition from the magnetized phase, c Free energies. Top: 
near the first-order phase transition at ^i/Ep = 1 for 0,2/Ep = -0.35, -0.1 and 0.15 for the black, 
blue and red traces respectively, as marked by the red flags in b. Bottom: near the second-order 
phase transition at ^2/Er = -2.5 for Vti/Ep = 4.5,5.5, and 6.5 for the black, blue and red traces 
respectively. 
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Figure 2 | Measured phase transition a Tensor magnetization Mzz measured as a function 
of Oi, showing both second-order [02(^^i = 0) = -2.0000(3)£^r] and first-order [02(^^i = 0) = 
-I.OE'r] phase transitions in comparison with theory. These curves followed the nominally hori¬ 
zontal trajectories (see Methods Summary) marked by red dashed curves in c. b Tensor magne¬ 
tization measured as a function of O2 at VIi/Er = 1.86(6), 1.48(4), 1.01(3), and 0.41(1), plotted 
along with the predicted critical ^2. In a and b the light-colored region reflects the uncertainty in 
theory resulting from our rs 5% systematic uncertainty in Oi. c Phase transition. Black (red) sym¬ 
bols depict data obtained using vertical (nominally horizontal) cuts through the phase diagram. 
Top, measured phase transitions plotted along with theory: solid (phase transition), and green 
vertical line (tricritical point, Of). Bottom: 20 % to 50 % transition width showing the clear shift from 
first to second-order with increasing Oi. 
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Figure 3 | Quenching dynamics The system was prepared in the unmagnetized phase with 
Oi = 0.74(8 )^h and ^^2 was ramped through the phase transition at ramp-rates dVt 2 /dt = 
-0.2, -0.3, -0.4, and -{).^Er/s (blue, black, red, and green symbols, respectively). The inset 
shows the decreasing width of the first-order transition as the ramp-rate decreases. 
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Figure 4 | Metastable states Top, Measured magnetization plotted along with theory. The sys¬ 
tem was prepared at the desired O 2 = -2Er, was then increased to its displayed final value; 
during this ramp Q 2 also changed, and the system followed the curved trajectory in the bottom 
panel. Each displayed data point is an average of up to 15 measurements, and the colored region 
reflects the uncertainty in theory resulting from our 5% systematic uncertainty in Oi. Bottom, 
state diagram: theory and experiment. Blue: two states; black: three states; white: one state. Col¬ 
ored areas denote calculated regions where the color-coded number of stable/metastable states 
are expected. Symbols are the outcome of experiment. Each displayed data point is an average 
of up to 20 measurements. 
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